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Abstract 

Systems with a first integral (i.e., constant of motion) or a Lyapunov function 
can be written as "linear-gradient systems" x = L(x)W (x) for an appropri- 
ate matrix function L, with a generalization to several integrals or Lyapunov 
functions. The discrete-time analogue, Ax/ At = LVV where V is a "discrete 

gradient," preserves V as an integral or Lyapunov function, respectively. 
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I. INTRODUCTION 



Integrals and Lyapunov functions — quantities that are conserved or dissipated, 
respectively — are fundamental in dynamics. They severely constrain the system's evolution 
and can be used to establish stability. There is no universal method to find such quantities, 
but if they are known (e.g., on physical grounds) we show that the system can be presented 
in a universal form which makes the conservation/dissipation property manifest. Although 
elementary, this result is very general and will find many applications: here we use it to 
preserve the conservation/dissipation property under time discretization. 

We start with the definition and an example of each of the classes of systems covered in 
this Letter. 

(i) Hamiltonian systems. 

Hamiltonian systems are ubiquitous in physics [|TJ. They have the form x = JW(x), 
x G M. 2n , where V(x) denotes the Hamiltonian function, and J := ( _j d ! q ), where Id 
denotes the identity matrix in IR n . 

Example 1. A simple Hamiltonian system is the pendulum |T[ x,\ = X2, %2 — ~ sin(xi); 
here n = 1 and V(xi,x 2 ) = \x\ — cos(xi). 

(ii) Poisson systems. 

Poisson systems also occur very frequently in phys-ics ( Q, App. 14). They have the 
form x = fl(x)'VV(x), x G M. n , where V(x) again denotes the Hamiltonian function 
and the Poisson structure Q(x) is an antisymmetric matrix (Q t (x) = —Q(x)), satisfying 
the Jacobi identity Qj^dk^igm + Qekdh^mj + ^mkdk^je. — 0. 

Example 2. The equations of motion of a free rigid body with moments of inertia 
I\, I2, and ^3 form a Poisson system with angular momentum x G M 3 and Poisson 
structure 
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Q{x) 



\ 



-£3 x\ 
x 2 —x\ 



(1) 



and Hamiltonian V(x) = \ Xlj=i x 1/h- (Actually this is an example of a so-called 
Lie-Poisson structure, in which Q(x) is a linear function.) 

(iii) Systems with a first integral. 

The ordinary differential equation (ODE) x = f(x), x G M n , is said to have the first 
integral V if dV{x)/dt = 0. 

Example 3. A Lotka-Volterra system ||. The ODE 



x l = e X3 , x 2 = e Xl + e X3 , x 3 = Be Xl + e X2 



(2) 



where B is a parameter, possesses the integral 



V(xx,x 2 , x 3 ) = e X2 Xl + B(x 2 - xi) - x 3 . 



(3) 



(iv) Gradient systems. 

Gradient systems arise, e.g., in dynamical systems theory |§. They are described by 
x = -VV{x), x e W 1 . 

Example 4- The system x\ = — 2x\(x\ — l)(2xi — 1), x 2 = —2x 2 is a gradient system 
[|], with n = 2 and V(x\, x 2 ) = x\{x\ — l) 2 + x\. 

(v) Systems with a Lyapunov function. 

The ODE x = f(x), x G M n , is said to possess the Lyapunov function V if dV(x)/dt < 
0. These functions were introduced by Lyapunov |13| and are a crucial ingredient of 
his direct or second method in the study of dynamical stability P,|l8 . 

Example 5. 
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Xl = —%2 — x\, X2 = X\ — x\ 



(4) 



has the Lyapunov function 



V(xi,x 2 ) = xj + x\. 



(5) 



What do the above five classes of dynamical systems have in common? A preliminary 
answer would be that they all possess a function V{x) such that dV(x)/dt < 0. That is, 
classes (i), (ii) and (iii) each possess a function V(x) such that dV(x)/dt = 0, and classes 
(iv) and (v) each possess a function V(x) s.t. dV(x)/dt < 0. 

In Section 2 we announce the result that classes of systems (i) to (v) have even more in 
common: under some mild technical assumptions, they can all be written as special cases of 
the novel class of "linear-gradient systems." In Section 3 we show how these linear-gradient 
systems can be integrated numerically in such a way that V(x) is constant or non-increasing, 
as appropriate. 

An extended version of this work, including proofs of the results presented here, is given 



in 14 . 



II. "LINEAR-GRADIENT SYSTEMS 



Our main result is the following: 



Theorem 1. Let the ODE 



x = f{x), 



f e C r , 



(6) 



possess a C r+1 Morse function V(x), where 



(a) 



dV 
~dt 







i.e., V is an integral; or 



(b) 




i.e., V is a (weak) Lyapunov function; or 



(c) 



dV 
~dt 



< where f\x) ^ 0, i.e. 



V is a strong Lyapunov function. 



4 



Then for all {x \ W(x) 7^ 0} there exists a locally bounded C r matrix L(x) such that the 
ODE (||) can be rewritten in the linear- gradient form 

x = L(x)VV(x), (7) 

where 

(a) L(x) is an antisymmetric matrix; resp. 

(b) L(x) is a negative semidefinite matrix; resp. 

(c) L(x) is a negative definite matrix. 
Remarks: 

1. A Morse function is a function whose critical points are all nondegenerate. A negative 
semidefinite matrix L is a matrix such that v f Lv < for all vectors v. A negative 
definite matrix L is a matrix such that v l Lv < for all non-zero vectors v. 

2. Under a coordinate transformation x \— > C(x) we have L(x) i— ► L(x) := 
dC(x)L(x)(dC(x)y. This implies that the theorem is invariant under coordinate trans- 
formations, because L is antisymmetric, negative semidefinite, resp. negative definite 
iff L is. 

3. The theorem has a converse: if an ODE is in linear-gradient form (0) with L antisym- 
metric, resp. negative semidefinite, resp. negative definite, then V is an integral, resp. 
weak Lyapunov function, resp. strong Lyapunov function. 

4. If the sign of dV/dt (zero, nonpositive, or negative) depends on x, then L can be 
chosen to be antisymmetric, negative semidefinite, or negative definite respectively, 
depending on x. The type of representation is not unique: at points where dV/dt = 0, 
L can be chosen to be either antisymmetric or negative semidefinite. 
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5. A particular L(x) satisfying the requirements of the theorem is 



Lij(x) 



fjVj ~ Vjfj + $ij fkVk 

E«2 



(8) 



where Vj = dV/dxj. However, L in (0) yielding @ is not unique. In particular, 
under further mild technical conditions, there is an L which extends smoothly through 
critical points of V. 

6. The fact that all systems with an integral can be written in the skew-gradient form 
x = L(x)VV(x) was, as far as we know, first published in [16]. The general case is 



new, although the special case of the converse with L(x) symmetric negative definite 
is well known and forms the subject of "generalized gradient systems" in dynamical 
systems 0. 



The constructive proof of Theorem 1 is given in [14| . We now give some illustrative examples 
of the above theorem. 

Example 6. Particle in ID with friction |J. Consider the ODE 

df{x l ) 



x x — x 2 , x 2 



dx\ 



- ax 2l 



(9) 



VV(x 1 ,x 2 ) 



(10) 



where a > is a coefficient of friction, and / is a potential function. Eq. (||) has the energy 
V(xi,x 2 ) = \x\ + f(xi) as a Lyapunov function, and can be written in the linear-gradient 
form (0) as 

x\ 1 

— 1 —a 

For a = 0, the system is conservative, and the matrix L is antisymmetric (case (a) above). 
For a > 0, the system is dissipative, V is a (weak) Lyapunov function, and L(x) is negative 
semi definite (case (b) above; cf §9.4 of 0). 

Example 7. An averaged system in wind-induced oscillation 0. Consider the system 



V 



/ 



X\ 



-C,x\ - \x 2 + x\_x 2 



x 2 = Xxi - (x 2 + ~{x\ 



. 'L- o 



(11) 



Here ( > is a damping factor and A is a detuning parameter. Guckenheimer and Holmes 
H remark that for ( = 0, Eq. (|TT| ) has the Hamiltonian 



V(xi,x 2 ) 



A 



1 



? (xl + xl) + - {^dxl - -± 



(12) 



and for A = 0, Eq. ( |TT| ) is a gradient system with 



z 2 ) = i (x? + a|) + 1 - 



;r 2 



(13) 



We now show that, for all allowed values of ( and A, Eq. (|TT|) can be written in linear- 
gradient form. To this end denote ( = pcos(6) and A = psin(#). Then Eq. ([LTD can be 
written in linear- gradient form i = LVV with 



^ —cos(6) — sin(#) ^ 



V 



sin(0) — cos(#) 



(14) 



and 



V(x u x 2 ) = (xl + xl) - ^ sin(0) faria| - y 



+ - COs(#) ( — x\x 2 \ ■ 



(15) 



Note that for the matrix L in this example we have v t Lv = — cos(#)|t> 2 |. Therefore, in the 
physical regime (where £ > and hence cos(#) > 0) either the matrix L is antisymmetric 
and V is an integral (for cos(#) = 0), or L is negative definite and V is a strong Lyapunov 
function (for cos(0) > 0). (Note that for A = £ = 0, we have p = and we are free to choose 
9. In this limit, the system possesses an integral 



Vi(x 1 ,x 2 ) 



as well clS cl Lyapunov function 



V 2 (x 1 ,x 2 ) 



1 / 2 X x 

2V 2 3 



1 /" C 2 2 

2ly- XiX2 



(16) 



(17) 



and V given by Eq. (|T^) represents an arbitrary linear combination of these two functions). 



Example 8. Reproduced from [[Uj, here is the linear-gradient form for the ODE 
Example 3: 

\ 



in 



X 2 

K x 3/ 



—er 



V 



w, 





-e Xl - 

e X3 e Xl _|_ gX 3 

where V is given by Eq. ([$]). 

Example 9. Here is the linear-gradient form for the ODE (H) in Example 5: 



(18) 



/ 



/ • ^ 




( 










a 


V 2 ) 




\ 


—b a ^ 



VV{xi,x 2 ) 



(19) 



where a = — {x\ + x^)/{x\ + x%), b = —(x\ + x\ + x%x\ — xix^)j{x\ + a;|), and V{xi,x-i) is 
given by Eq. (|5]). Note that the matrix L in (|l^) is negative definite. 



III. DISCRETE GRADIENTS AND THE NUMERICAL INTEGRATION OF 

LINEAR-GRADIENT SYSTEMS 



For differential equations whose time evolution has particular structural properties, such 
as preservation of symplectic structure, phase space volume, symmetries, or conserved quan- 



tities, it is desirable to mimic these properties in any numerical integration | 20| . This is 
particularly useful in long-time integrations. One can also view the discrete-time analogues 



as interesting physical systems in their own right |L2 . 

A major application of the linear-gradient formulation (|7|) is that it has a simple and 
elegant discrete-time analogue; moreover, this analogue is also a universal representation for 
systems of each class. 

Definition 1. Let V(x) be a differentiable function. Then W(x, x') is a discrete 
gradient of V if it is continuous and 



W(x, x') ■ (x' -x) = V(x') - V(x), 
VV(x,x) = W(x). 
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(20) 



Discrete gradients are not unique. Several examples of discrete gradients are given in 

0B0- 



Definition 2. The function V is an integral of the map x t— > x' if V(x') = V(x), Vx. 
It is a weak Lyapunov function if V(x') < V(x), Vx. It is a strong Lyapunov function if 
V(x') < V(x) for all x such that x ^ x' . 

Theorem 2. [JlJ] Let the map x i— > x' be defined implicitly by 

=) = L{x, x', r)VV(x, x') (21) 

where VV is any discrete gradient, L is a matrix function, and r represents a timestep. Then 
V(x) is an integral, resp. weak Lyapunov function, resp. strong Lyapunov function of the 
map if L is antisymmetric, resp. negative semidefinite, resp. negative definite. Conversely, 
for any map with such a V, and any discrete gradient V, at points such that W(x, x') ^ 
there exists an L such that that map takes the form ft2l\). 

It follows that is a discrete approximation to the linear-gradient system (0) that 
preserves integrals, resp. Lyapunov functions, provided the method is consistent, i.e., 
L(x, x, 0) = L(x). 

Equations similar to fl2Hp, ( |2~T| ) have appeared in many energy-conserving schemes for 
Hamiltonian systems [0,f|,[l(],|lll] , although the first axiomatic presentation was and the 



first application to all systems with an integral was |17 |. 



IV. CONCLUDING REMARKS 

(i) In this Letter, for simplicity we have restricted our discussion to the case of one first 
integral or Lyapunov function. In |14[ we show that an n-dimensional ODE with 



m < n — 1 integrals and/or Lyapunov functions Vi, . . . ,V m , can be written in the 
"multilinear-gradient" form 

x = L(x)VV l . . . W m , x G M n , (22) 
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where L(x) is an (m + l)-tensor. Structure-preserving integrators for Eq. (|2~2 ) have 
also been constructed, generalizing Eq. (|2~T|). 

(ii) Associated with (multi) linear-gradient systems of the form (p2] ) there is also a formu- 
lation in terms of a bracket: 

df(x) 



dt 

where the bracket is defined by 



U,V u ...,V m } L , (23) 



ii,...,t p 1 p 

where p = m + 1. This bracket satisfies the Leibnitz rule in each of its variables: 

{fl, fj-h 0(#1> ■ ■ ■ ,9k), fj+l, • • • , 

= (f~) {A'"--' /i-i' • • • ' ( 25 ) 

j = 1, . . . ,p. Conversely, the tensor L is defined by the fundamental brackets L ilt ^ jip = 
{x iiy . . . , Xi }l. It follows that V(x) is an integral, resp. weak Lyapunov function, resp. 
strong Lyapunov function of the (multi) linear-gradient system ([22]) iff W = Vx, resp. 
W < Var, resp. W < for all x such that | W(x)| + 0, where := {V, V u ... , V m } L . 
It also follows that V is an integral (resp. Lyapunov function) of the system (B^) iff 
Vj is an integral (resp. Lyapunov function) of the system 

x = L(x)VVi . . . VVj-iVVVVj+i . . . W m , (26) 

where 

L — L- ■ ■ (27) 



Special cases of the bracket fl2J) are the Poisson bracket and the Nambu bracket pi 



(iii) We hope to address the numerical order of accuracy of the integrator Q2"T| ) in a forth- 
coming publication. 
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